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Abstract 

We consider scalar perturbations in the time-dependent Horava- 
Witten Model in order to probe its stability. We show that during 
the non-singular epoque the model evolves without instabilities until 
it encounters the curvature singularity where a big crunch is supposed 
to occur. We compute the frequencies of the scalar field oscillation 
during the stable period and show how the oscillations can be used to 
prove the presence of such a singularity. 
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1 Introduction 



The classical Big-Bang model based on Einstein's General Relativity does not 
explain the origin of the initial singularity, which is expected to be understood 
at a higher dimensional string theory model. 

It is known that the Horava-Witten theory [1] relates 11-dimensional su- 
pergravity on the orbifold S 1 /Z 2 with strongly coupled heterotic E 8 x E 8 
string theory. Then, this link suggests that as one probes higher energies, 
our 4-dimensional universe would first go through an intermediate regime 
where the orbifold dimension becomes visible, thus making our world appear 
5-dimensional [2, 3], and only at energies of order the string scale the universe 
would look 11-dimensional. 

The intermediate 5-dimensional energy regime has led to the conjecture 
that our universe is a brane in a higher dimensional world. These Horava- 
Witten inspired "braneworlds" have a topology defined by a line times a 
non-compact space with two branes at the boundaries of the line element, 
which is usually characterized as the orbifold S 1 /Z 2 . Due to the compactness 
of the orbifold direction, one brane has positive tension while the other has 
the opposite negative one. By construction these spacetimes are supersym- 
metric and this unbroken supersymmetry of the static background solution 
guarantees the local stability of the model [4] . 

However, the situation changes when the background solution becomes 
time-dependent. The idea of having moving branes has naturally led to 
the suggestion that the hot big bang is the result of a brane collision [5], 
or even that our universe undergoes an endless sequence of cosmic epochs 
which begins with a universe expanding from a big bang and ending with a 
contraction to a big crunch [6]. 

Colliding braneworld scenarios offer an alternative point of view for many 
long-standing problems and puzzles such as those connected with the exis- 
tence of gravitational shortcuts (apparently faster than light signals) [7]. 
These issues refer to the different perspective of higher dimensional theories 
and conventional 4-dimensional ones. In fact, the proposed observational 
tests of braneworld scenarios must be able to distinguish between higher 
and lower dimensional effects. Moreover, the former effects must appear as 
corrections of the latter in these models. 

One of the approaches used to study these collisions is through low energy 
supergravity limits. Since we lack a complete M-theory formulation, the 
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solutions come from the 5-dimensional supergravity theory obtained after 
a dimensional reduction of the higher dimensional scenario. In the case of 
Type II-B theory exact solutions to this problem have been found in [8], 
where the motion of the D3 branes towards one another leads to the complete 
disappearance of the universe in a spacetime singularity. 

In the context of heterotic brane models an exact solution of the equations 
of motion derived from the same Horava-Witten lagrangian has been found 
[9]. As in the previous case, the setup obtained from a dimensional reduction 
to 5 dimensions evolves to a singularity where the entire spacetime is annihi- 
lated. This can be taken as an indication of an inherent classical instability 
in braneworld models consisting of a positive and a negative tension brane. 
However, it has recently been proposed a way to avoid the annihilation by 
introducing an arbitrarily small amount of matter on the negative tension 
brane [10]. This alternative scenario is specially appealing for the ekpyrotic 
and cyclic universe models [5, 6] due to the momentary disappearance of the 
orbifold dimension caused by the collision and the subsequent bounce after 
which both branes continue to expand into the future. 

In this paper we study the stability of the solution in [9] by introducing 
a scalar perturbation. General perturbations of solutions to the Einstein 
equations are very well known in the literature, see e.g. [11] and further 
references therein. In the present case we deal with perturbations of time- 
dependent solutions, a more delicate issue [12]. Although we know from 
previous discussions that there exists an instability [8, 9], our motivation is 
to check the possibility of obtaining some information about the signatures 
we would expect before the spacetime collapses. 

The paper is organized as follows. In section 2 we describe the time- 
dependent Horava-Witten solution under study. In the next section we intro- 
duce the scalar field perturbation, we solve the corresponding Klein-Gordon 
equation in the bulk, and we discuss the field evolution in the setup. Finally, 
we present our conclusions in section 4. 

2 The Time-Dependent Horava-Witten Model 

One of the most popular models connecting M-theory with phenomenology is 
the Horava-Witten Model [1]. In this setup the 11-dimensional spacetime is 
the product of a Calabi-Yau compact manifold and a 5-dimensional spacetime 
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consisting of two parallel 3-branes or domain walls, one with negative tension 
and the other with positive one. When we perform a dimensional reduction 
from eleven to five dimensions, we obtain a 5-dimensional supergravity theory 
which admits an exact static supersymmetric solution of the form [2] 

dsl = H(-dt 2 + dx 2 ) + H 4 dy 2 , (1) 

where 

H=l + k\y\, = -31ogtf, (2) 

and k 2 = 2m 2 / 3 is a constant, while is a scalar field characterizing the size 
of the internal Calabi-Yau space. The equations of motion for the metric and 
can be obtained from the Lagrangian 

C b = ^-g(R- l -(d<P) 2 -m 2 e 2 <t>). (3) 

In this scenario a second domain wall is introduced at y = L, taking y as 
a periodical coordinate with period 2L and making the Z 2 identification, 

y <-»• -y- 

The model based upon this solution is not physically realistic because it 
is static, while our universe is time-dependent. Led by this motivation many 
attempts to introduce dynamics in this model have been made. Remarkable 
examples are the Ekpyrotic Model [5] , where the big bang is regarded as the 
result of the collision of an external brane with our universe, and the Cyclic 
Universe [6], where the distance between the branes oscillates in time. 

However, recently Chen et al. [9] found that the lagrangian (3) also 
admits a time-dependent 3-brane solution given by 

ds 2 5 = H l l 2 (-dt 2 + dx 2 ) + Hdy 2 , (4) 

with 

H = ht + k\y\ , = --log#, (5) 

where k 2 = 8m 2 / 3 and h is an arbitrary constant. 

If we turn off the time-dependence in this solution, i.e., set t — 0, we 
can go back to the time-independent solution (l)-(2) by performing the co- 
ordinate transformation y = |y 2 and rescaling the constant Co of the original 
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work of Lukas et al. [2], where H = k\y\ + cq. In fact, this transition be- 
tween the solutions can be better visualized if we set Co = 0, which is always 
possible since it is an integration constant. 

The solution (4)- (5) represents a bulk spacetime with two 3-branes with 
negative and positive tensions as boundaries at y = and y = L, respec- 
tively. The interval < y < L is of S 1 /Z 2 type, where S 1 occupies the strip 
—L < y < L and the Z 2 symmetry corresponds to the identification y <-> y. 
By choosing h < as in [13] we notice that H is a positive but decreasing 
quantity when t < 0. The spacetime is bounded in the range < y < L, 
but the proper length of this interval is time-dependent. Thus, the separa- 
tion between the branes decreases as t becomes positive, i.e., the universe is 
contracting. When t turns out to be positive, H vanishes along —ht = ky. 
This line represents a singularity which first develops on the negative tension 
brane (y — 0) at t — and travels along the orbifold dimension reaching 
the positive tension brane (y = L) after a finite time t = kL/(—h). The 
spacetime can not be extended beyond this point. Thus, before the branes 
effectively collide the power-law singularity wraps the whole spacetime. 

Although we know that an instability exists in the present model, we still 
want to check the possibility of describing its existence before the spacetime 
collapses. 



3 The Scalar Perturbation 

We consider a perturbation in the form of a scalar field obeying the massive 
Klein-Gordon equation 

□ $ = m 2 $ ; (6) 

or equivalently 

1 d M (^g MN d N )^ = m 2 ^. (7) 



Using the metric (4) in spherical coordinates and the derivatives of H, 
d t H = h and d y H = fcsgn(y), Eq.(7) can be rewritten as 



X 



1 „ , . s 1 



sin v sin 



d e (smed ) + — —Jl + H~ l d 2 y + -H-hgn{y)d y - m 2 U> = . (8) 
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Now we decompose the scalar field $ as 

${t,r,9,(f),y) = Z(t,r,y)Y £m (6,4>) , 
where the spherical harmonics obey the equation 



sm sin 9 v 

Thus, the function Z(t,r,y) satisfies the equation 

- Aaz + d 2 z + -d r z - e -^±^-z + h-w&z + 



H 



+ -H- 3/2 sgn(y)d y Z - m 2 H 1/2 Z = 




100 



100 -loo 



Figure 1: Potential a 2 + m 2 H 1 ! 2 for a 2 = 1, m — 0.1. 
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In order to separate the radial part we write Z as Z(t, r, y) = \l/(t, y)R(r). 
Substituting back into (11) we obtain 



d 2 r R + \r + [a 2 - ^J^j R = 0, (12) 

and 

+ _ -l= d ^ _ ^f^* + {a 2 + m 2 i^ 2 )* = o , (13) 

where a 2 is a constant. 

The solution to equation (12) is given by 

Rir) = jF J (l + e - a ') + 7f Y (l + e - ar )- (14) 

where J and Y are Bessel functions of first and second kind, respectively, 
and A and B are integration constants. 

The solution of equation (13) must be obtained by numerical methods. 
We must thus introduce some values for the parameters h, k, a 2 , and m. 
Following [13] we choose h to be negative so that our description of the 
system begins at a negative time and evolves to a positive time. In order to 
simplify the calculation we choose the simplest positive values for k and a 2 , 
k = \h\ and a 2 = 1, as well as a positive scalar field mass m. 

The potential term (a 2 + m 2 H 1 / 2 ) is shown in Fig.(l). We see that it 
is positive-definite although it possesses some singular derivatives along the 
line \y\ = t. 

Our results for the evolution of the field \I/ are shown in Figs. 2-4. 

Fig. 2 provides information on the stability of the system before t — 0, 
when the singularity develops on the negative tension brane. We can see that 
for earlier times the system is stable. The oscillations show a quasinormal 
phase where we can calculate quasinormal frequencies (see Table 1), which 
show a strong dependence on the value of the parameters \h\ and k when 
m^O. We notice that the real part of the frequency increases with \h\ and 
k, as well as with m. Since we calculated the frequencies in the stability range 
where the oscillations are almost uniform, the imaginary part appears to be 
very small, and its sign just suggests a slight change in the amplitude of the 
oscillation (positive sign indicates an increase, while negative sign indicates 
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Figure 2: Evolution of the massless (left) and massive (right) scalar field \I/ 
at the position of the negative tension brane y = for several values of h. 



Table 1: Quasinormal Frequencies at y = 0. 



|h| 


m = 


m = 0.1 


m = 2 


k 














1 


1.428 


-0.0023 


3.452 


0.0015 


6.411 


0.0010 


2 


1.428 


-0.0021 


4.028 


0.0016 


7.662 


0.0025 


3 


1.428 


-0.0021 


4.363 


0.0014 


8.491 


0.0026 


4 


1.428 


-0.0020 


4.689 


0.0017 


8.976 


0.0022 


5 


1.428 


-0.0020 


4.909 


0.0016 


9.520 


0.0028 



a decrease). However, as soon as the scalar field detects the formation of the 
singularity, the oscillations increase in amplitude and frequency, revealing 
the instability of the system. 

Once the singularity has formed on the negative tension brane, it trav- 
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Figure 3: Evolution of the massless (left) and massive (right) scalar field \I/ 
at a bulk position y = L/2 between the two branes for several values of h. 



Table 2: Quasinormal Frequencies at y = L/2 = 50. 



|h| 


m = 


m = 0.1 


m = 2 


k 














1 


1.293 


-0.0005 


3.740 


0.0006 


7.140 


0.0013 


2 


1.288 


-0.0005 


4.363 


0.0008 


8.491 


0.0010 


3 


1.293 


-0.0006 


4.760 


0.0008 


9.240 


0.0017 


4 


1.293 


-0.0006 


5.150 


0.0003 


10.134 


-0.0025 


5 


1.293 


-0.0006 


5.417 


0.0007 


10.472 


0.0022 



els along the extra dimension in the direction of the positive tension brane. 
Analizing our results at a point between the two branes as shown in Fig. 3 
we observe that the frequency of the scalar field oscillation increases grad- 
ually until encountering the singularity. Table 2 shows some quasinormal 
frequencies calculated by taking the most stable oscillations. 

After a time t = kL/(—h) the singularity reaches the positive tension 
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Figure 4: Evolution of the massless (left) and massive (right) scalar field ^ 
at the position of the positive tension brane y = L for several values of h. 



Table 3: Quasinormal Frequencies at y = L = 100. 



|h| 


m = 


m = 0.1 


m = 2 


k 






^R 




^R 




1 


1.199 


0.0076 


3.927 


0.0773 


7.854 


0.0863 


2 


1.200 


0.0639 


4.620 


0.0642 


9.240 


0.0799 


3 


1.204 


0.0582 


5.150 


0.0590 


10.472 


0.0838 


4 


1.213 


0.0548 


5.512 


0.0551 


10.833 


0.0232 


5 


1.213 


0.0519 


5.818 


0.0524 


11.220 


0.0271 



brane. In Fig. 4 we see how unstable the system becomes at this point. The 
quasinormal frequencies shown in Table 3 provide evidence of this behaviour, 
since their positive imaginary part becomes larger than in the previous cases. 
The singularity eventually wraps both branes and the space between them, 
which constitutes the bulk, smearing the actual brane collision. 
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4 Conclusions 



In this paper we have considered the time-dependent Horava-Witten solution 
found by Chen et al. [9] . We studied the stability of the model by introducing 
a massive scalar perturbation. 

We found that before any singularity appears, i.e., for t < 0, the model 
evolves without any instability. For such a time-dependent solution it was 
possible to calculate quasinormal frequencies at different positions along the 
extra coordinate. We chose 3 representative points to show our results, 
namely, at the position of each brane and between the two branes. The 
quasinormal frequencies show a strong dependence on the parameters of the 
model for a massive scalar field. The larger the parameters \h\ and k are, or 
alternatively the larger the mass of the scalar field m is, the larger the real 
part of the fundamental frequency is. This dependence disappears when the 
mass of the scalar field vanishes. 

At t — a curvature singularity forms on the negative tension brane. This 
singularity travels along the extra dimension and hits the positive tension 
brane at t — kL/(—h), before the two branes actually collide. 

From t — to t — kL/(—h) the scalar field oscillations increase frequency 
and amplitude showing the instability generated by the curvature singularity 
that finally envelopes all the spacetime. 

Despite knowing the instability of this model that leads to a big crunch, 
we have shown that it is possible to detect it. The evolution of the scalar field 
used for this purpose can provide some information about the parameters of 
the model and the extension of the extra dimension. 

A complete picture of the instability should also include gravitational 
perturbations. However, since both the metric and the potential are time- 
dependent, such a procedure is beyond the aims of the present paper. 
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